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ABSTRACT 

Experiments  are  described  that  reveal  the  three-dimensional 
character  of  the  plastic  zone  in  front  of  notches  and  cracks  in 
plates  of  an  Fe-3Si  steel  and  a  plain  carbon  steel.  These  define 
the  plane-stress  regime  as  a  function  of  applied  stress  and  plate 
thickness.  They  also  provide  a  rationale  for  the  DM  (Dugdale- 
Muskhelishvili)  model  as  a  tentative  elastic-plastic  solution  of 
a  crack  under  plane  stress.  Refinements  that  offer  a  way  of 
taking  work  hardening  and  rate-sensitive  plastic  deformation  into 
account  are  described.  In  this  way,  unnotched  tensile  properties 
—  the  stress-strain  curve  a-;d  reduction  in  area — are  used  to 
calculate  plastic-zone  size,  crack-tip  displacements  and  strains, 
the  crack-extension  stress,  and  the  fracture  toughness,  in  accord 
with  experiments.  Finally,  the  approach  is  extended  to  ductile 
crack  propagation  and  used  to  calculate  the  crack  speed  and  the 
stress,  strain,  and  strain  rates  imposed  on  material  in  advance 


of  a  moving  crack. 
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INTRODUCTION 


(1) 

Treatments  of  crack  extension  in  metals  based  on  the  solutions  of  Inglis 

(2) 

and  Griffith  fail  because  they  regard  the  metal  as  completely  elastic.  The 

modified-elastic  solutions  of  Orowan^and  Irwin^are  more  successful.  These 

make  adjustments  for  local  plastic  deformation,  but  they  do  not  describe  stress 

and  strain  in  the  region  adjacent  to  the  crack  tip  Further  progress  demands 

elastic-plastic  solutions  that  can  describe  in  detail  what  happens  in  the 

plastic  zone,  and  these  will  make  it  possible  to  predict  crack  extension  from 

the  unnotched  mechanical  propjrties.  Inroads  into  the  elastic-plastic  problem 

(5)  (6-9) 

have  been  made  by  Allen  and  Southwell,  and  others,  and  further  progrees  can 
be  expected.  Yet,  at  this  writing,  there  are  no  rigorous  solutions  that  define 
the  local  stress  and  strain  for  a  crack  in  tension  under  either  plane  strain 
or  plane  stress;  the  solution  for  the  intermediate  case  is  not  even  in  sight. 

This  paper  explores  two  alternate  methods  of  attacking  the  plane-stress 
problem:  (1)  experimental  measurements  and  (2)  a  special  elastic-plastic 

solution  which  compromises  generality  and  rigor  for  simplicity.  The  experi¬ 
ments  reveal  the  character  of  the  plastic  zone  in  front  of  cracks  and  notches 

in  an  Fe-3Si  steel  and  a  plain  carbon  steel  plate.  In  this  way,  they  offer 

(10,11) 

a  rationale  for  the  DM  (Dugdale-Muskhelishvili)  model.  It  has  already  been 

(10,12) 

shown  that  this  model  makes  reasonably  good  predictions  of  plastic-zone 
size  and  crack-tip  displacement  under  plane-stress  conditions.  In  this  paper, 
a  refined  version  of  the  DM  model,  capable  of  taking  work  hardening  and  rate- 
sensitive  yielding  into  account,  is  described.  Plastic-zone  size,  crack-tip 
displacements  and  strains,  the  crack-extension  stress,  and  Kc^  are  calculated 
in  this  way  from  unnotched  properties.  The  calculations  are  also  extended  to 
the  case  of  a  propagating  ductile  crack,  and  they  define  the  crack  speed  as 


t  Linear  elastic  fracture-toughness  parameter. 


2 

well  as  the  stresses,  (trains,  and  strain  rates  generated  in  advance  of  the 
crack.  Where  possible,  calculations  are  compared  with  experiment,  and  the 
resulting  correlations  support  the  usefulness  of  the  DM  model. 


OBSERVATIONS  OF  PLASTIC  ZONES 


Procedure  and  Materials 

This  section  deals  with  the  three-dimensional  character  of  the  plastic 

zone  in  front  of  a  crack  and  how  it  changes  with  nominal  stress  and  plate 

(12) 

thickness.  The  experiments  follow  along  lines  previously  reported,  but  are 

more  complete  and  systematic.  The  plastic  zones  were  generated  by  edge  notches 

in  rectangular  coupons  of  an  Fe-3Si  steel  (Si,  3.317.;  C,  0.047.)  and  a  plain 

(13) 

carbon  steel  (Project  Steel  E:  C,  0.227.;  Mn,  0.367.).  These  coupons,  derived 
from  stress-relieved  plates  (1  hour  at  475°C)  of  warm-rolled  Fe-3Si  steel^ 
and  hot-rolled  Steel  E,  were  machined  in  thicknesses  ranging  from  0.453  inch 
to  0.017  inch. 

Two  types  of  notches  were  studied:  (1)  a  0.25-inch-deep  by  0.006- inch¬ 
wide  saw  cut  (root  radius  ~0. 003  inch),  and  (2)  a  0.25- inch-deep  fatigue 
ttt 

crack.  After  machining  and  notching  (and  fatiguing),  the  coupons  were 

tttt 

recrystallized  to  anneal  out  traces  of  deformation  introduced  during  machining 

and  fatiguing.  One  set  of  Fe-3Si  steel  coupons  was  annealed  at  1200°C,  a  high 

temperature  intended  to  coarsen  the  grain  size  and  eliminate  discontinuous 

yielding.  As  a  final  step,  the  surfaces  of  the  test  coupons  were  electro- 

(14) 

polished  with  the  Morris  reagent. 


t  Eight  inches  long  by  2.5  inch  wide,  with  a  centrally  located 
4- inch-long  gage  section. 

tt  Reduced  50  percent  at  300°C. 

t+t  Grown  from  0.12-inch-long  by  0.006- inch- wide  saw  cuts  by 

cycling  the  coupons  in  tension  between  4,000  and  38,000  psi. 

tttt  Fe-3Si:  800°C  for  1  hour;  Steel  E:  900°C  for  1  hour. 
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(b>  Static  Yield  Stress 

(c)  For  this  steel.  Equation  B-2  was  replaced  by 


4 


The  notched  coupons  were  then  slowly  loaded  to  various  stress  levels  in 
a  tensile-testing  machine;  the  peak  load  was  maintained  for  4  minutes  and  then 
gradually  released.  Stress  levels  are  reported  either  as  T,  the  net  section 
stress,  or  as  T/Y,  the  ratio  of  net  section  stress  to  yield  stress.  As  the 
next  step,  residual  strains  normal  to  the  plate  surface  were  detected  and 
recorded  with  the  aid  of  an  interference  microscope.  Following  this,  Fe-3Si 
steel  samples  were  aged  20  minutes  at  150°C  to  decorate  the  dislocations. 

Plastically  deformed  regions  on  the  surface  of  Fe-3Si  steel  were 
revealed  by  etch-pitting  electrolytically  in  Morris  reagent.  The  coupons 
were  then  sectioned,  repolished,  and  etched  to  reveal  the  plastic  zone  in 
the  interior.  The  plain  carbon  steel  could  not  be  etched.  In  its  case, 
plastic  zones  are  revealed  on  the  surface  by  virtue  of  slight  surface  tilts 
which  accompany  the  plastic  deformation. 

True-stress/true-strain  curves  of  the  Fe-3Si  and  plain  carbon  steels 
are  reproduced  in  Figure  la,  and  other  data  are  given  in  Table  1.  Figure  1 
shows  graphically  the  stress-strain  formulations,  described  in  Table  1,  that 
serve  as  inputs  to  the  computer  calculations  and  are  discussed  later.  A 
complete  list  of  symbols  and  their  definitions  appears  in  Appendix  A. 

Experimental  Results 

Figure  2  shows  the  plastic  zones  revealed  on  the  surface  of  Fe-3Si  steel 

plates  of  different  thicknesses  and  their  development  with  increasing  applied 

(12) 

stress.  The  interpretation  given  previously  is  that  there  are  two  different 
types  of  zones — one  associated  with  plane  strain,  and  the  other  with  plane- 
stress  conditions: 

1,  Plane  Strain.  At  low  stress  levels,  two  yielded  regions  fan 
out  in  directions  roughly  normal  to  the  plane  of  the  crack  (the 
±y-direction,  Figure  2g;  see  Figures  2c,  2f,  2g,  and  2h);  at  this 


7  (True  Stress ) ,  1000  psi 


6 


Fig.  2.  Influence  of  Nominal  Stress  and  Plate  Thickness  on 

the  Plastic  Zone  Revealed  by  Etching  the  Surfaces  of 
Machine-Notched  Fe-SSi  Steel  Coupons ( 800° C  Anneal) : 

(a)(b)  0.017 -inch  thick,  c/t=l4.  7 

(c)(d)(e)  0.058-inch  thick,  c/t=  4.4 

(f)(g)(h)&(i)  0.200-inch  thick,  c/t=  1.3 

Morris  reagent,  oblique  lighting,  magnification  67.  5X. 


<»)T/r»0.4  (b)  T/Y  ■  0. 6  > 


(0  T/Y  *0.4  <g)  T/Y  *0.6  (h)T/r»0.8  (1)  T/Y -0.8 


fans.  On  the  surface,  the  two  wedges  are  always  separated  by  a 
distance  corresponding  to  the  plate  thickness  (see  Figures  2b,  2e, 
and  2i).  Etching  interior  sections  parallel  to  the  plate  surface 
reveals  that  the  two  wedges  merge  into  a  single  wedge  on  the  mid- 
section,  as  shown  in  Figure  3.  Normal  sections  (see  Figure  4) 
confirm  that  the  wedges  are  in  reality  wedgelike  regions  inclined 
«t  45  degrees  to  the  tensile  axis  vdiere  shear  has  occurred  through 
the  thickness  direction.  At  this  stage  there  is,  in  fact, 
noticeable  strain  through  the  thickness  direction  in  the  form  of 
a  depression  (dimple)  on  the  surface. 

More  detailed  studies,  now  in  progress,  reveal  45-degree-shear  near  the  plate 
surface  and  in  the  interior,  even  when  the  plastic  zone  is  much  smaller  than 
the  plate  thickness.  It  is,  therefore,  more  correct  to  regard  Fe-3Si-steel 


Fig,  4.  Appearance  of  Plastic  Zone  on  the  Front  Surface ,  Back  Surface ,  and  a 
Normal  Section 

(Fe-ZSi,  machine  notched ,  800° C  anneal)  0.1 9? -inch  thick,  T/Y  =  0.0) 
Morris  reagent,  oblique  lighting ,  magnification  8.ZX, 

plastic  zones  as  composites:  partly  hinge,  and  partly  inclined  wedge.  At 
low  stress  levels,  particularly  for  short  cracks  and  thick  plates,  the  hinge 
component  predominates,  while  for  suitable  combinations  of  high  stress,  long 


crack,  and  thin  section,  the  inclined  wedges  make  the  major  contribution. 


The  plastic  zones  observed  on  the  electropolished  surfaces  of  the  plain 
carbon  steel  (Figure  5)  are  also  interpreted  as  inclined  wedges.  Here, 
although  a  number  of  wedge-shaped  zones  are  nucleated,  the  growth  of  only  one 
is  favored,  and  this  wedge  is  longer  than  the  corresponding  wedges  in  the 
Fe-3Si  steel  (compare  Figures  2  and  5  and  note  the  difference  in  magnifica¬ 
tion)  .  Since  the  wedges  appear  above  the  crack  centerline  on  one  side  of  the 
plate,  and  below  on  the  opposite  surface,  they  undoubtedly  represent  a  single 
45-degree-inclined  wedge. 

It  can  be  seen  in  Figures  6a  and  6b  that  notch  acuity  has  little  effect 

on  either  size  or  the  shape  of  the  inclined  wedge  component,  that  is,  the 

tt 

zone  character  under  plane-stress  conditions.  Also  shown  in  Figures  6c  and 
6d  (these  correspond  to  Curves  1.1  and  1.2  in  Figure  la)  is  the  effect  of  the 

“if  There  can  probably  be  interpreted  as  Luders'  bands. 

ft  There  is  a  marked  effect  of  notch  acuity  on  the  distribution  of 
strain  within  the  hinge  component,  particularly  close  to  the 
root  radius.  This  will  be  described  in  more  detail  in  a  future 
publication. 


I 


Fig.  6.  Plastic  Zones  Produced  by  Loading  0.060-inch- thick  Coupons  to  T/Y=  0.78 
and  Revealed  on  the  Surface  by  Etching  and  Interferometry : 

(a)  Fe-ZSi,  machine  notched ,  800° C  anneal  (Curve  1.1 ,  Fig.l ),  etched 

(b)  Fe-3Si,  fatigue  cracked ,  800° C  anneal  (Curve  1.1,  Fig.l),  etched 

(c)  Same  as  (a),  interferometry 

(d)  Fe-3Si,  machine  notched,  1200° anneal  (Curve  1.2,  Fig.  1),  etched 
(c)  Steel  E,  machine  notched  (Curve  1.3,  Fig.  1),,  interferometry 

(a)  &  (b)  Oblique  lighting.  In(c)  &  (d),  percent- .strain  values  are 
quoted.  In(e) ,  the  numbers  give  relative  displacements  in  units 
of  0.00001  inch.  The  closely  spaced  contours  -2,  0,  +2,  +6  are 
consistent  with  a  45-degree-inclined  wedge-shaped  region  sheared 
about  1%. 


Magnification  12, 2X. 


yield  drop.  Comparison  reveals  that  the  difference  between  a  material  with  a 

yield  drop  and  one  without  is  less  than  experimental  error.  This  evidence 

.  „  (16)  (17) 

refutes  the  suggestion  of  Durelli  et  al,  and  Dixon  and  Visser  that  the 


wedge-shaped  plastic  zone  is  due  to  the  yield  drop,  as  such.  It  also  opens 

(18) 

to  question  the  view  o£  McClintock  and  Iryin  that  the  effect  is  mainly  the 
result  cf  a  low  work-hardening  rate.  As  a  contrast  to  the  Fe-3Si  steel  contour, 
Figure  6e  shows  the  profile  for  plain  carbon  steel.  Comparison  with  the 
.visual  observations  (for  example.  Figure  3)  reveals  that  the  zone  represents 
a  tilt  on  the  surface,  as  suggested  previously. 

The  cransition  from  a  dominant  hinge  to  a  dominant  wedge  (for  example, 
from  plane  strain  to  plane  stress)  is  gradual  and  therefore  difficult  to 
define.  It  seems  likely  that  the  hinge  will  dominate  so  long  as  45-degree- 
shear  is  constrained;  that  is,  so  long  as  the  inclined  wedges  are  contained  by 
elastic  material =and  do  not  penetrate  to  the.  plate  surfaces ..  Such  .penetration 
will  not  have  pccurred  so  long  as  the  portion  of  the  hinge  visible  on  the 
surface  extends  a  distance  pjj  <t/2  on  either  side  of  the  crack  (see  Figure  2g; 
t  is  the  plate  thickness) .  This  serves  as  an  approximate  criterion  for  plane 
strain.  The  contribution  of  the  wedges  probably  begins  to  exceed  that  of  the 
hinge  when  the  wedge  length  p  >4t  (see  Figure  2i) .  This  is  merely  a  rough 
estimate  and  serves  as  a  tentative  criterion  for  plane  stress.  In  both 
instances,  pH  and  p  increase  with  increasing  applied  stress.  Consequently, 
the  transition  from  plane  strain  to  plane  stress  will  occur  at  higher  stresses 
in  thicker  sections.  This  is  shown  qualitatively  in  Figures  2  and  5  and  more 
quantitatively  in  Figure  7,  For  example,  the  0.017-inch-thick  coupon  (Figure 
2a)  already  displays  45-degree-shear  prominently  when  the  nominal  stress  is 
40  percent  of  the  yield  stress;  whereas  in  the  0.200- inch- thick  coupon, 
stresses  in  excess  of  90  percent  of  yield  are  required  (Figure  2i) . 

From  the  evidence  presented,  it  is  possible  to  construct  stylized  three- 
dimensional  pictures  of  the  plastic  zone  under  plane-stress  conditions.  As 
shown  in  Figure  8,  the  zones  displayed  by  the  plain  carbon  and  Fe-3Si  steels 
are  basically  similar.  They  consist  of  either  a  long,  single,  45-degree- 
inclined  wedge,  or  two  long  45-degree- inclined  wedges  intersecting  at  90 
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Fig.  9.  The  DM  Model:  (a)  Crack  with  Wedge-shaped  Plastic  Zones,  (b)  Crack  with 

" Equivalent "  Internal  Tension ,  (c)  Uniform  Tension  Distribution  of  Dugdale, 
(d)  Varying  Tension  Distribution,  and  (e)  Stepwise  Approximation  of  the 
Varying  Tension  Distribution. 


degrees.  Figures  8d  and  8c  show,  schematically,  that  the  wedges  are  attached 
to  the  hinge  portion;  but  the  exact  pattern  of  deformation  in  the  region  where 

the  hinge  and  wedge  merge  is  not  yet  clear. 

THE  DM  MODEL 


The  Simple  Model 

In  the  light  of  the  preceding  section,  Dugdale' a  model  of  a  crack  with 
wedge-shaped  plastic  zones  takes  on  special  significance.  The  model,  illus¬ 
trated  in  Figures  9a  and  9b,  consists  of  a  crack  of  length  2c  in  an  infinite 


plate.  When  the  plate  la  subjected  to  a  nominal  stress,  I,  plastic  zones  of 
length  p  are  generated  at  either  end.  Together,  the  zones  and  the  crack  are 
assumed  to  have  the  shape  of  the  stressed  elastic  slit  shown  in  Figure  9b. 

The  material  within  the  plastic  zone  is  then  replaced  by  an  "equivalent" 
distribution  of  internal  tension,  S(x),  acting  on  the  slit  interface — equiv¬ 
alent  in  the  sense  that  it  simulates  the  support  that  was  derived  from  the 
yielded  material.  This  final  step  transforms  the  elastic-plastic  problem 
into  a  purely  elastic  one  that  can  be  treated  by  the  Muskheliahvili  conformal 

.  ’  (H)t 

mapping  procedures. 

The  DM  model  possesses  a  zone  lying  in  the  plane  of  the  crack,  as 
opposed  to  the  45-degree- inclined  wedge(s)  that  dominate  real  plane-stress 
-zones  (see  -Figure  8c)-i  The  model  also:  possesses-  a  special  wedge-shape,, 
effectively  a  Tresca  yield  criterion,  and  a  material  that  undergoes  no  fur¬ 
ther  elastic  deformation  after  yielding.  These  assumptions  may  not  be  overly 
restrictive.  For  example,  the  distinction  between  a  straight  and  an  inclined 
wedge  is  minimized  as  p  becomes  larger  than  t,  and  vanishes  as  t  -*  0.  The  DM 
model  is  useful — to  the  extent  the  various  assumptions  are  applicable  or  not 
crippling — because  it  transforms  the  complex  three-dimensional  plastic  region 
into  a  much  simpler  one-dimensional  zone. 

,  (1°) 

Dugdale  considered  the  simple  case  vrtiere  S(x)  is  a  uniform  distribution 
equivalent  to  Y,  the  yield  stress  of  the  material  (see  Figure  9c).  For  this 

case,  the  following  can  be  formulated  explicitly  in  terms  of  c,  T,  and  Y: 

(12) 

(a)  The  elastic  stress  gradient,  o(x)  for  y  ■  0 

(10) 

(b)  The  plastic-zone  size,-  p 

(20) 

(c)  The  displacement  gradient,  v(x)  for  c  <  x  £  p 

(20) 

(d)  The  crack-tip  displacement,  vc  =  V(x  «  c) . 

The  relevant  equations  are  summarized  in  Table  2, 


t  Barenblatt  'has  adopted  a  similar  approach  to  treat  the  finite 
cohesive  strength  of  completely  brittle  materials. 


TABLE  2.  IMPORTANT  EQUATIONS  OF  THE  DM  MODEL 


5 •SZ ?  J'i  L'V.l  >-»-• 
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So  far,  albeit- limited  experimental  tests  of  these  equations  have  been 

(10) 

encouraging.  Dugdale's  zone-size  measurements  performed  on  mild  steel  are 

(12) 

described  very  well  by  Equation  4.  Hahn  and  Rosenfield  found  that  zone  size 
and  crack-tip  displacement  values  measured  for  Fe-3Si  steel  are  in  reasonable 
accord  with  Equations  4  and  5a.  For  example,  the  criterion  for  plane  stress, 
p  >  4t,  can  be  formulated  via  Equation  4  as 

-1 


t  >* 


sec 


TT  T 

2Y 


1 


(14) 


and  this  is  consistent  with  the  experimental  measurement  given  in  Figure  7 . 
The  criterion  for  plane  strain,  pjj  <  -r,  can  be  formulated  on  the  basis  of 

(21)  i 

Tetelman  s  observation  that  pjj  ■j,  and  is 


i  < 


n  T  ' 
sec  Yy  -  1 


-1 


(15) 


and  this  is  also  in  reasonable  accord  with  the  measurement.  Equations  14  and 
15  indicate  that  the  zone  character  is  determined  by  the  dimensionless  param¬ 
eter,  c/t,  and  this  given  the  results  of  Figure  7,  obtained  solely  by  changing 
t,  greater  generality. 

At  the  same  time,  some  discrepancies  between  theory  and  experiment  have 

(10,12) 

been  noted.  At  the  higher  stresses,  Dugdale's  measured  p  values  and  those 

for  Fe-3Si  steel  (see  Figure  7  for  c/t  <  2,  where  some  of  the  half-filled 
points  lie  in  the  plane-stress  region)  are  consistently  smaller  than  predicted. 

This  is  related  to  work  hardening;  the  increase  with  strain  of  o,  the  flow 

stress  in  the  plastic  zone,  means  that  the  value  of  S  becomes  larger  than  Y. 

At  the  same  time,  the  load-bearing  cross  section  within  the  plastic  zone 

diminishes  as  a  consequence  of  constant-volume  plastic  deformation,  and  this 

ultimately  leads  to  local  necking  and  a  reduction  in  S.  The  net  effect  is  an 
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Modified  DM  Model 

Work  hardening  and  constanC-volume  deformation  influence  two  things: 

(1)  the  distribution  of  internal  tension,  S(x),  and  (2)  the  y  distribution  of 
strain  in  the  plastic  zone,  which  is  shown  in  Figure  B-l  and  discussed  more 
fully  in  Appendix  B.  This  distribution  fixes  the  relation  between  v  and  7, 


7  =  f(v) 


(16) 


where  v  is  the  displacement  at  a  given  distance,  x,  along  the  abscissa,  c  £  x 
fi  (c  +  p),  and  7  is  the  maximum  true-plastic  strain  at  the  same  distance.  At 
present,  the  form  of  Equation  16  must  be  derived  from  experiment,  and  the 
simple  model  presented  in  Appendix  B  offers  an  approximate  way  of  formulating 
the  distribution  in  terms  of  the  measurable  quantity,  d,  the  width  of  the 
plastic  zone  (see  Figure  5f), 


(16a) 


To  simplify  matters  further,  the  ratio  d/t  is  assumed  to  remain  constant. 

The  idea  of  relating  strain  and  displacement  in  terms  of  a  fixed  distance,  d, 

(20) 

has  been  advanced  by  Goodier  and  Field. 

The  same  model  also  describes  t",  the  minimum- section  thickness  at  a 
given  distance,  x,  along  the  abscissa,  c  £  x  £(c  +  p),  and  this  connects  the 
value  of  S  with 

s  =  c  (“ )  =  a  exp  (-7)  ,  (17) 

and  (via  Appendix  B), 


S 


d  -  2v 
d  +  2v 


(17a) 


The  flow  strength  and  work-hardening  rate  of  material  within  the  zone  is 
expressed  by  the  true-stress/true-strain  characteristics, 


a  -  f(r) 


(18) 
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Some  o£  Che  calculations  Chat  will  be  described  exploit  the  following  a  -  IT 
equation, 

o  -  Jnie+hf  ,  (18a) 

where  Y,  B,  and  n  are  material  constants.  This  is  equivalent  at  all  but  small 

n 

strains  to  the  more  usual  power  law,  "5  =  KT  ,  but  has  the  advantage  that 
C  (e  =  0)  =  Y . 

Together,  Equations  5,  16,  17,  and  18,  or  the  more  specific  forms  used 
here  (Equations  16a,  17a,  and  18a),  offer  a  first  estimate  of  internal-stress 
distribution,  S(x) .  The  resulting  function  is  approximated  by  a  stepwise 
distribution,  as  shown  in  Figure  9e.  Then,  via  Equations  7  to  13  in  Table  2 
(these  are  the  analogies  to  Equations  1  to  6  for  the  case  of  a  stepwise 
internal-stress  distribution)  and  particularly  Equation  12,  successive  ap¬ 
proximations  of  S(x)  can  be  obtained.  This  reiterative  process  leads  to  the 
correct  form  of  S(x)  and  appropriate  values  of  p,  <j(x),  v(x),  and  vc,  as  well 
as  ¥(x)  and  e^.  The  procedure,  while  prohibitive  by  hand,  is  easily  managed 
on  a  digital  computer.  The  flow  diagram  of  a  computer  program  is  given  in 
Appendix  C.  It  is  programmed  for  an  IBM-7094  Fortran  II  system  located  at 
Wright-Patterson  Air  Force  Base,  Ohio. 

Plastic-Zone  Calculations 

Results  of  computer  calculations  exploiting  the  modified  DM  model  and 
results  derived  from  the  simpler  form  are  presented  in  Figures  10  and  11. 

The  values  in  Figure  10  were  calculated  for  Curve  1.4  (Figure  lb,  Table  1) 
intended  to  match  a  4330M  steel.  Figure  10a  compares  the  distribution  of 
internal  tension;  Figure  10b  shows  the  elastic  stress  gradient  and  true- 
stress  values  within  the  plastic  zone;  Figure  10c  gives  corresponding  values 
of  the  plastic  strain  gradient.  It  should  be  noted  that  the  strain  values 
quoted  at  the  various  relative  distances,  x/c,  are  the  true  strains  at  the 
minimum  section  (see  Equation  B-2,  Appendix  B);  the  value  of  the  crack-tip 
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strain,  e-^  s  =TC  (see  Figure  10c),  is  the  maximum  strain  generated 

within  the  DM  zone.  Figures  10b  and  10c  illustrate  that  the  crack  not  only 

concentrates  the  stress,  but  is  a  potent  "strain  concentrator"  as  well. 

Figure  11a  shows  that  the  plane-stress  zone  size  measurements  of  Dugdale 

(p/t  ^  4),  while  in  accord  with  the  simple  theory,  agree  even  better  with 

the  predictions  of  the  modified  DM  model.  Figure  lib  shows  similar  results 

for  Fe-3Si  steel.  This  is  good  evidence  that  the  systematic  deviation  from 

the  simple  model  can  be  attributed  to  work  hardening  rather  than  a  strain- 

(10) 

rate  effect.  The  correlation  also  lends  some  confidence  to  the  other 
calculations  which  have  not  been  tested. 

Figures  12a  and  12b  suranarize  plastic-zone  calculations  for  the  hypo¬ 
thetical  "high-strength"  and  "conventional"  steels  in  Figure  lc  (Curves  1.7, 
1.8,  and  1.9;  Table  1).  These  are  intended  to  show,  in  a  general  but  quanti¬ 
tative  way,  the  influence  of  strength  level  and  work-hardening  rate  on  zone 
size  and  crack-tip  strain.  Figure  12a  illustrates  that  the  strength  level 

has  little  effect  on  zone  size  when  the  comparison  is  made  at  the  same 

+ 

relative  nominal-stress  level,  T/Y,  This  follows  from  Equation  4a,  valid  at 
low  stress.  On  the  other  hand.  Figure  12b  shows  that  strength  level  can  have 
a  significant  effect  on  Tc, even  at  comparable  T/Y  values;  the  crack  produces 
larger  strains  within  the  higher  strength  material*.  This  has  important  impli¬ 
cations  for  fracture  and  is  a  consequence,  first,  of  the  factor  Y/E  which 
appears  in  Equation  5a  and,  secondly,  of  Equation  B-2,  for  example, 


Thus,  the  crack-tip  strain  also  depends  on  plate  thickness.  The  value  of  the 

work-hardening  exponent  has  a  modest  effect  on  p  and  ?c>  via  its  effect  on 

the  S(x)  distribution.  The  quantity  n  may  actually  have  a  larger  influence 

than  is  at  first  apparent,  because  it  probably  modifies  the  zone  width,  d. 
t  Assuming  n  and  d/t  are  unchanged, 

tt  At  constant  Y  and  d/t. 
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Fig.  12.  Influence  of 
strength  level , 
work  hardening, 
end  ductility  as 
shown  by  calcu¬ 
lations  for  the 
hypothetical 
steels  of  Fig.  1c 
(Curves  1.7,  1.8, 
and  1.9,  and 
Table  1) 
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APPLICATION  OF  THE  MODIFIED  DM  MODEL 
TO  CRACK  EXTENSION 


In  the  preceding  section,  the  DM  formalism  was  used  to  treat  the  growth 
of  a  wedgelike  plastic  zone.  The  model  can  also  deal  with  the  growth  of  the 
crack  itself  if  the  cracking  is  by  ductile  shear  or  fibrous  fracture.  These 
mechanisms  usually  involve  large  strains  locally,  and  the  plane-stress 
criterion,  p  >  4t,  can  be  satisfied.  The  analysis  can  take  advantage  of  the 
fact  that  states  of  stress  and  strain  generated  within  the  wedgelike  DM  zone 
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are  similar  to  those  produced  in  the  neck  of  an  unnotched  plate  or  sheet 
t 

coupon.  Indeed,  fractures  in  unnotched  tensile  bars  of  copper  involve 

(22) 

thin  zones  of  heavily  deformed  material.  The  cracks  form  within  the  zones  at 

(22) 

a  terminal  stage  of  hole  coalescence  ,  irrespective  of  vrtiether  the  plate  is 

notched  or  unnotched,  and  this  stage  is  conveniently  identified  by  the  strain 

ft 

at  fracture.  In  this  way,  a  critical  strain  criterion  for  fracture  of 
material  in  the  region  adjacent  to  the  crack  tip  (the  most  heavily  strained 
region)  can  be  adopted, 


ec  ^ 


_* 


(20) 


_* 


where  e-  is  the  ductility  of  the  unnotched  material  in  terms  of  true  strain 
as  calculated,  for  example  from  RA,  the  reduction  in  area  [  T  =  -/n(l-RA)]. 
Equation  20  is  also  the  criterion  for  crack  instability  (crack  extension), 


because  the  strain  imposed  on  material  ahead  of  the  extending  crack  continues 

d  rc 

dc 

* 


to  increase,  >  0  (Equations  5a  and  B-2) ,  while  its  capacity  for  strain, 


T  ,  remains  constant.  The  quantity  T*  h  is  t^ie  nominal  crack 

extension  or  fracture  stress  of  the  cracked  plate,  and  Kc  *  c  »  where 
Kc  is  the  linear  elastic  fracture- toughness  parameter.  Values  of  T*  and  Kc 
can  be  computed  with  the  same  equations  and  program  described  in  the  last 
section.  For  the  given  conditions; 

(a)  Geometry,  c  and  t 

(b)  Flow  properties,  $  =  f(T)  (for  example,  Y,  B,  and  n,  Table  1) 
and  d 

jl 

(c)  Ductility,  e  (reduction  in  area), 


+  The  large  triaxial  stress  component  generated  when  the  crack  is  under 
plane  strain  is  absent  under  plane  stress  (by  definition).  The  smaller 
triaxial  stress  component  associated  with  the  neck  is  present  in  both 
cases.  The  conditions  in  the  notched  coupon  are  also  similar  to  those 
in  a  round  test  bar,  since  the  influence  of  the  biaxial  stress  component 
(normal  to  tensile  axis  and  in  the  plane  of  the  plate)  on  the  stress  and 
strain  at  fracture  is  probably  small. 

t+  Provided  the  rates  of  straining  are  comparable,  the  fracture  condition 
is  represented  by  a  single  point  on  the  a  -  s’  curve.  Either  the  true 
stress  or  the  true  strain  at  fracture  identifies  this  point;  but 
strain  is  more  convenient. 
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and  an  assumed  value  of  T,  the  quantity  ec  is  calculated  and  compared  with 


T  is  then  appropriately  adjusted  and  the  procedure  repeated  until 
fc  -*  r  ;  concurrently,  T  -*  T*.  Two  exceptions  to  this  procedure  are  noted 


in  the  footnote  below. 

Calculations  of  this  kind  were  tailored  for  a  4330M  steel  and  a  2219-T87 


aluminum  alloy.  The  crack-extension  stress  of  these  materials  had  previously 

(23) 

been  measured  over  a  range  of  crack  lengths  (see  Figure  13) .  Fractures  were 


of  the  ductile  shear  variety.  The  a  -  e  curves  could  be  constructed  from 


published  data,  and  RA  values  are  reported  in  the  literature  (see  Figure  lb. 


Curves  1.4  and  1.5,  and  Table  1).  Only  the  values  of  d,  the  zone  width,  have 


not  been  measured;  the  values  quoted  in  Table  1  are  rough  estimates  based  on 


observations  of  similar  materials. 


As  shown  in  Figure  13,  both  the  general  trend  and  the  absolute  values  of 
T"  and  Kc  calculated  in  this  way  are  in  good  accord  with  the  experimental 


values.  It  is  particularly  significant  that  the  modified  DM  theory  predicts 

the  decrease  in  Kc  evident  when  y  >  0.8,  an  effect  t.iat  is  beyond  the  reach 

(23) 

of  linear  elastic- fracture  mechanics.  Figure  13a  shows  curves  calculated 


for  two  widely  differing  values  of  d/t,  and  this  gives  some  indication  of 


the  uncertainty  introduced  by  incomplete  knowledge  of  this  parameter. 


The  influence  of  strength  level,  work-hardening  rate,  and  ductility  on 


T*  and  Kc  was  calculated  for  the  two  hypothetical  steels  of  Figure  lc,  and 

t  Exception  1:  Failure  by  yielding  at  T  =  Y.  In  some  cases,  ec  <  ^  even 
when  Ti  Y,  Failure  is  then  assumed  to  occur  by  gross  yielding  and 
T*  =  Y,  An  alternative  approach,  not  used  here,  is  to  consider  the 
plate  at  a  uniform  strain,  ^T  =  ^ (g  =  T) •  ^  crack  of  length  c  is  now 

cut,  and  the  plastic  zone  is  the  region  where  s’  >  ex.  The  origin  of 
the  a  -  s’  curve  is  shifted  so  that  ex  corresponds  to  s  =0,  and  Y  = 
o r—  —  \.  The  critical  strain  is  nows'  =  s’  -  st>  and  the  model  can 

be  applied  as  before. 

Exception  2:  Failure  by  yielding  at  T  <  Y,  In  some  instances,  ec  <  e  , 
but  the  ratio  p/c  becomes  very  large.  The  failure  is  then  a  form  of 
plastic  instability,  and  T*  2  Tp  _ 


ft  For  constant,  d/t. 
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Fig.  IS.  Comparison  of  Calculated  Fracture  Stress  and  Ka  Levels  with  Values 
Measured  experimentallu . 


a.  Absolute  Effects  b.  Relative  Effects 


Fig.  14.  Calculated  Curves  Showing  the  Influence  of  Held-Strength  Level, 

Work-Hardening  Hate ,  and  Ductility  on  the  Fracture  Stress  of  Hypothetical 
Steels  (Equation  1 ,  Table  1). 

the  results  are  summarized  in  Figures  12c  and  12d,  All  three  material  param¬ 
eters  play  a  significant  role  in  determining  fracture  strength.  Figure  12c 
illustrates  that  there  is  not  a  one-to-one  correlation  between  strength  level 
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and  the  fracture  stress;  the  relative  fracture  stress  T*/Y  at  30,000  psi  is 
higher  than  at  the  300,000-psi  strength  level  (see  Figure  12c).  This  effect, 
which  assumes  importance  for  all  "high-strength"  materials,  has  as  its  origin 
the  role  of  Y/E  in  Equation  5b.  The  net  result,  shown  more  clearly  in  Figure 
14a,  is  that  the  AT1f  associated  with  a  fixed  AY  decreases  as  the  strength 
level  is  raised  when  T”,  n,  and  d  are  constant.  In  practice,  some  duc¬ 
tility  and  work-hardening  capacity  are  usually  sacrificed  to  obtain  higher 
strengths,  with  the  result  (as  shown  by  the  shaded  band  in  Figure  14a)  that 
the  fracture  stress  becomes  relatively  insensitive  to  yield  strengths  above 
200,000  psi.  This  means  that  the  benefits  (in  terms  of  increased  fracture 
toughness)  of  strength-level  increases  above  200,000  psi  are  likely  to  be 
marginal  if  they  are  accompanied  by  decreased  ductility  and  work-hardening 
rate. 

The  kind  of  "trade-offs"  that  may  be  encountered  for  a  300,000-psi  steel 
(see  Curve  1.7,  Figure  lc,  and  Table  1)  can  be  anticipated  by  the  calculations 
of  the  type  summarized  in  Figure  14b,  For  example,  a  33  percent  or  100,000- 
psi  increase  in  Y,  accompanied  by  a  20  percent  loss  in  both  the  RA  and  nt 
would  net  only  a  modest  9  percent  or  17,000-psi  gain  in  fracture  stress. 
Calculations  of  this  kind  offer  useful  guidelines  to  the  metallurgical  origins 
of  fracture  toughness  and  the  potential  value  of  specific  material  improve¬ 
ments  . 

APPLICATION  OF  THE  MODIFIED  DM  MODEL  TO 
CRACK  PROPAGATION 

The  balance  of  forces  between  the  elastic  matrix  and  the  plastic  zone 
for  stresses  below  T'v  is  a  form  of  static  equilibrium.  When  the  stress 
exceeds  T*  and  the  crack  begins  to  accelerate,  the  system  enters  a  state  of 
dynamic  equilibrium  controlled  by  the  dynamic  yield  characteristics,  provided 

t  This  corresponds  to  a  change  from  RA  =  30%  (e"  =  0.35)  to 
RA  =  24%  (e"  =  0.27),  and  from  n  =  0.05  to  n  =  0.04. 
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the  plastic  response  of  the  material  is  rate  sensitive^  The  main  difference 
is  that  the  distribution  of  internal  tension  now  reflects  the  flow  stress 
consistent  with  both  the  local  strain  and  the  local  strain  rate.  In  this 

*n» 

way,  the  rate  with  which  the  material  strains  to  the  critical  value,  e"  , 
determines  the  crack  speed.  When  the  effective  crack  length  grows  continu¬ 
ously,  the  equilibrium  is  constantly  displaced  and  the  crack  accelerates. 

In  special  cases,  the  effective  crack  length  is  limited  by  a  finite  dimension 
of  the  structure,  and  a  steady  state  may  be  attained.  This  mechanism  does 

not  consider  inertia  effects  and  is  only  valid  for  crack  speeds  and  acceler- 

(24) 

ations  substantially  below  the  limits  set  by  the  dynamic  analysis  of  Mott 

(25)tt 

and  Roberts  and  Wells. 

When  the  inertia  effect  can  be  ignored,  the  problem  can  be  treated  with 

ttt 

the  modified  DM  model.  The  only  additional  requirement  is  a  description  of 
the  true  stress  in  terms  of  strain  and  strain  rate, 

a  =  ffi,  ?)  .  (21) 

A  specific  example  of  Equation  21,  derived  from  dislocation-dynamics  consider- 

(26) 

ations,  is  uspful  for  this  purpose: 

.1  .1 

o  =  Ax  en  +  Ce)m  [Ag  +  A3  e]  m  .  (22) 

It  can  describe  an  upper  and  a  lower  yield  point  (see  Curves  1.10,  Figure  Id), 
different  degrees  of  rate  sensitivity,  and  its  constants  can  have  theoretical 
significance.  The  equation  is  used  to  perform  the  steady-state  calculations 


t 

tt 


ttt 


If  the  resistance  to  plastic  deformation  increases  with  increasing 
strain  rate.  (24)  (25) 

The  velocity  calculated  by  Mott  and  Roberts  and  Wells  >  relevant 
for  the  rate- insensitive  material. 


u  =  0.38  us  Jl  -  , 

can  be  regarded  as  an  upper  limit  (us  is  the  longitudinal  wave 
velocity,  and  c*  is  the  crack  length  at  instability). 

Goodier  and  Fielcl  (lave  applied  the  DM  model  to  the  dynamic  problem, 
but  their  calculations  do  not  take  rate-sensitive  yielding  into 
account. 
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described  in  the  following  paragraphs.  These  are  of  a  preliminary  nature, 
and  are  intended  merely  to  demonstrate  the  application  of  the  DM  model  to  the 
dynamic  case. 

The  crack  speed  and  the  stress,  strain,  and  strain-rate  gradients  in 

front  of  the  crack  were  calculated  for  a  propagating  ductile- shear  fracture. 

The  case  considered  is  a  crack  traveling  along  the  x-axis  at  steady-state  and 
t  * 

at  constant  velocity  for  a  fixed  nominal  stress,  T  >  T  ,  under  the  following 
conditions: 


(i)  Geometry .  ~  =  200 

(ii)  Dynamic  flow  characteristics.  Relation  III,  Table  1; 

the  constants  were  selected  to  approximate  the  response 

-1  +3 

of  mild  steel  in  the  strain-rate  range  10  to  10  per 

second  (see  Curves  1.10,  Figure  Id,  and  Table  1).  There 

is  evidence  that  the  constants  used  slightly  overestimate 

rate  sensitivity  of  steel  when  t  <  10  ^  sec  ^  and  seriously 

x  3  -1 

underestimate  if  when  e  >  10  sec  .  Relation  III 
(Table  1)  is  reasonable  for  a  mild  steel  with  Y  = 

_4  _i 

40,000  psi  when  measured  at  e  =10  sec  ,  and  this 
must  be  taken  into  account  in  future  calculations. 

The  expression  for  d  (footnote  in  Table  1)  is  derived 
from  actual  measurements  of  a  shear  fracture  in  a  mild 
steel  involving  a  fast-moving  crack. 


(iii)  Ductility .  The  value  e  =  0.33  (RA  =  287.)  was  derived 
from  the  shear  fracture  mentioned  above;  the  value 
e*  =  0.40  (RA  =  337.)  was  for  purposes  of  comparison. 

The  numerical  procedure  begins  with  a  trial  stress  distribution  and  a 

(1  6  * 

trial  velocity, u0.  First,  approximations  of  the  quantities  e,  e,  S  , 

t  Essentially  a  moving  axis  analysis,  since  the  effective  crack  length 
must  remain  constant. 
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and  S(x)  are  obtained^  and  the  process  is  repeated  until  the  solution  con- 
verges.  Then  e"c  is  compared  with  e  ,  the  trial  value  uD  is  adjusted  accord- 

M  4MV  JU 

ingly,  and  the  process  of  evaluating  ec  is  repeated.  In  this  way,  as  ec  -*  e  , 
u0  -*  u,  and  stress,  strain,  and  strain-rate  gradients  consistent  with  the 
velocity  are  determined.  The  singularity  of  strain  gradient  at  x  =  a  is 
avoided  in  the  computation  by  smoothing  out  the  strain-distribution  curve 
near  the  elastic-plastic  boundary.  This  process  is  justified  because  the 
cusp  is  very  localized  in  nature,  and  in  reality  one  does  not  expect  such  a 
cusp  to  occur. 

Figures  15  and  16  are  examples  of  the  results  of  dynamic  calculations. 

Figure  15  shows  the  relation  between  nominal  stress  and  crack  velocity.  The 

velocity  appears  to  be  exceedingly  sensitive  to  T  for  speeds  below  500  feet 

3  -1 

per  second.  Higher  speeds  generate  strain  rates  in  excess  of  10  sec  in 

the  plastic  zone  (see  Figure  16c),  and  the  steel  probably  enters  the  regime 

of  greater  rate  sensitivity.  For  this  reason,  the  crack  velocity  is  likely 

to  be  less  sensitive,  and  this  is  shown  qualitatively  by  the  dashed  portions 

of  the  curves  in  Figure  15.  The  calculations  also  indicate  that  ductility 

can  have  a  significant  influence  on  crack  speed. 

Stress,  strain,  and  strain-rate  gradients  are  reproduced  for  two  crack 

speeds  in  Figure  16.  The  stress  gradients  reveal  peaks  at  the  leading  edges 
of  the  plastic  zone  which  are  manifestations  of  the  yield  point.  These  are 
in  doubt,  and  are  shown  in  dashed  lines  for  two  reasons:  (1)  the  computa¬ 
tions  are  inaccurate  in  this  region,  and  (2)  there  is  some  question  whether 

the  deficiency  of  mobile  dislocations,  that  accounts  for  the  yield  point  in 

(26) 

an  annealed  material,  can  exist  at  an  elastic-plastic  interface.  The 

600-ft/sec  velocity  was  chosen  as  an  example  because  shear  fractures  travel- 

(27) 

ing  at  speeds  from  500-800  feet  per  second  are  encountered  in  practice. 

+  4  ~  d  T  d  e  _  . 

T  ®  =  uo>  o  is  given  by  Equation  22. 


c/t  *  200 
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Figure  16c  indicates  that  the  behavior  of  such  fast-moving  fractures  is 

controlled,  at  least  in  part,  by  the  material's  response  to  strain  rates 
2  4-1 

from  10  to  10  sec  .  This  range  is  relatively  unexplored — a  factor  limit¬ 
ing  crack-speed  calculations  of  this  kind.  The  high  rates  also  mean  that 
present  calculations  probably  underestimate  the  stresses  in  front  of  the 
600-ft/sec  crack  (Figure  16a) . 

The  present  calculations  tend  to  confirm  the  results  obtained  earlier 
(28) 

by  Hahn,  Reid,  and  Gilbert  with  a  much  cruder  model:  (1)  stresses  of  the 
order  of  100,000  psi  are  attained  ahead  of  fast-moving  cracks  in  steel,  and 
(2)  the  plastic  zone  is  smaller  at  higher  speeds  (Figure  16a) .  The  high 
stresses  plus  the  triaxial  component  associated  with  a  shorter  and  less- 
wedgelike  zone  will  favor  cleavage  fracture.  It  is  anticipated  that  further 
calculations  of  this  kind  will  prove  useful  for  describing  the  ductile-to- 
brittle  transition  behavior  of  steel. 


DISCUSSION 

The  DM  analysis  is  versatile,  and  there  is  evidence  supporting  its 
relevance  in  specific  cases.  The  real  questions  are  its  general  applica¬ 
bility  and  the  precision  this  approach  offers — and  these  are  points  that 
warrant  discussion  and  more  experimentation.  Clearly,  real  plastic  zones 
deviate  from  the  ideal  DM  wedge.  How  seriously  do  these  and  other  assump¬ 
tions  affect  the  predictions  of  the  model? 

Long  wedgelike  zones  have  now  been  observed  under  plane-stress  condi¬ 
tions  (Equation  14)  in  silicon  steel,  in  plain  carbon  steei^’^  "^ »  in 
(32)  (29,33) 

high-strength  steel  ,  and  in  cold-rolled  copper.  The  wide  range  of  con¬ 
figurations  employed  in  these  investigations  tends  to  confirm  our  observa¬ 
tion  that  the  appearance  of  the  DM  zone  does  not  depend  critically  on  notch 
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geometry.  Consequently,  this  type  of  zone  is  not  an  isolated  phenomenon. 

It  seems  likely  that  all  of  these  are  wedges  inclined  at  45  degrees  to  the 
tensile  axis.  This  deviation  does  not  appear  to  be  serious  in  view  of  the 
agreements  between  theory  and  experiment. 

In  some  cases,  wedges  also  intersect  the  plate  surface  at  45  degrees 

to  the  x-axis  (refer  to  Figure  9) .  These  "side  bands"  have  been  observed  by 

(17)  (16)  (33) 

Dixon  and  Visser  and  Durelli  et  al  in  steel,  and  by  Druyvesteyn  et  al  in 

hard- rolled  copper  sheet,  but  not  by  the  present  investigators.  In  general, 

they  are  found  in  conjunction  with  the  more  usual  zones  discussed  above. 

The  type  of  deformation  contributed  by  side  bands  and  their  influence  on  the 

local  strain  distribution  is  not  clear. 

(30) 

The  plastic  zones  observed  by  Bateman  et  al  on  the  surfaces  of  cracked 
aluminum  alloys  are  not  elongated  wedges  but  circular  in  character.  At  first 
glance,  they  may  seem  a  departure  from  the  conceptual  picture  developed  here; 
but  this  is  not  the  case.  Bateman's  tests  involve  combinations  of  stress, 
crack  length,  and  thickness  that  place  the  the  zone  in  the  lower  portion  of  the 
transition  region  between  plane  stress  and  plane  strain  (Figure  7) .  Contour 
maps  of  the  circular  zones  are  comparable  to  the  one  obtained  for  the  zone 
shown  in  Figure  2h.  This  does  not  prove  that  all  plane- stress  zones  are 
wedgelike,  but  again  shows  that  both  the  stress  level  and  the  crack  length- to- 
thickness  ratio  must  be  taken  into  account  in  comparisons. 

The  hinge  component  close  to  the  tip  of  the  crack  represents  another 
deviation  from  the  theoretical  DM  zone.  There  is  seme  evidence  that  the  be¬ 
havior  of  this  portion  of  the  plastic  zone  is  similar  to  the  wedge.  Bateman's 
(30)  (22) 

measurements  and  those  of  Tetelman  indicate  py,  the  extent  of  predominantly 

hinge  zones,  obeys  the  same  functional  relation  as  p  and  is  about  one-half  the 

value  predicted  for  plane  stress.  Furthermore,  the  equations  describing  zone 

(10)  (9) 

size  and  displacements  for  plane-stress  tension  and  both  plane-  and  anti-plane- 
(19) 

shear  are  nearly  identical.  In  spite  of  this,  the  hinge  complicates  crack 
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extension.  Within  the  influence  of  the  hinge,  stable  crack  growth  ("pop- in") 

(34) 

is  possible  along  lines  formulated  by  McClintock.  This  comes  about  because 
the  deformed  regions  extend  out  to  large  y-distances.  As  the  crack  extends, 
these  regions  are  left  behind  and  new  zones  take  their  place  at  the  crack  tip. 
A  crack  with  this  history  of  local  yielding,  extension,  and  more  local  yield¬ 
ing  (as  opposed  to  a  virgin  crack  of  the  same  length)  generates  smaller  strains. 
6  €  c  ^l. 

Consequently  g""  *  can  be  negative,  and  the  condition  ec  >  ”  may  not  be 

maintained.  Stable  crack  growth  within  the  hinge  in  the  early  stages  of 

loading  can  be  incidental  to  the  unstable  growth  mechanism  fostered  by  the 

wedge-shaped  zones;  but  it  will  modify  the  stress-strain  distribution,  and 

may  cause  departures  from  the  DM  model.  The  plastic-zone  history  will  also 

tend  to  stabilize  the  shear  crack  growing  within  the  wedge  zone.  In  this  case, 

plastic  deformation  is  largely  confined  to  the  region  immediately  in  front  of 

the  crack;  the  effect  is  thus  likely  to  be  smaller,  and  the  condition  ^-4-^  >  0 

Oc 

may  still  be  fulfilled.  Shear  crack  growth  implies  that  the  crack  is  in  the 
plane  of  a  45-degree-inclined  wedge,  and  this  presents  a  problem  since  the 
starting  crack  is  usually  pictured  normal  to  the  tension  axis.  The  crack  must 
therefore  reorient  itself,  and  this  probably  occurs  in  the  region  influenced 
by  the  hinge  component.  The  existence  of  this  reorientation  stage  need  not 
radically  alter  the  criterion  for  crack  extension  employed  in  this  paper,  but 
probably  modifies  it  in  some  way.  For  these  reasons,  agreements  cited  here 
between  crack-extension  experiments  and  the  DM  theory  should  be  regarded  with 
caution. 

Still  another  important  deviation  is  the  fact  that  the  zone  width,  d, 
and  displacement,  v,  do  not  coincide  (see  Figure  9a).  While  this  does  not 
invalidate  zone-size  or  zone-displacement  predictions,  it  does  preclude  the 
possibility  of  calculating  strain  directly.  This  is  a  crucial  point,  because 
the  distribution  of  strain  in  the  y-direction  is  a  prerequisite  for  calculat¬ 
ing  ductile  fracture  and  for  implementing  the  work-hardening  and  strain-rate 
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corrections.  The  concept  of  a  strain  distribution  involving  a  constant  width, 
d,  exploited  here,  is  convenient  but  oversimplified;  both  d  and  the  character 
of  the  distribution  may  change  with  stress  and  distance  from  the  crack  tip. 
Certainly,  a  general  solution  of  the  configuration  of  the  neck  in  an  unnotched 
sheet  or  plate  would  offer  useful  insights — but  this  problem  apparently  has 
not  been  worked  out.  Further  experiments  characterizing  the  strain  distribu¬ 
tion  under  load  may  offer  the  possibility  of  correlating  this  with  the  shape 
of  the  stress-strain  curve. 


CONCLUSIONS 

Under  plane-stress  conditions,  localized  plastic  zones  at  notches  and 
crack  in  an  Fe-3Si  steel  and  a  plain  carbon  steel  are  dominated  by  long  45- 
degree-inclined  wedge-shaped  regions.  These  wedges  resemble  the  plastic 

zones  attending  Dugdale's  plane-stress  model  (DM  model)  of  a  crack  in  an 
elastic-plastic  material.  The  mathematical  treatment  of  the  "simple"  DM 
model  offers  reasonable  descriptions  of  zone  sizes  and  displacements  in  real 
materials.  A  "modified"  form  of  the  same  model  can  deal  with  work-hardening 
and  rate-sensitive  plastic  deformation.  In  this  way,  the  stress,  strain, 
and  strain- rate  gradients,  the  crack-extension  stress,  effects  of  work  hard¬ 
ening  and  ductility,  and  the  speed  of  propagating  cracks  can  be  calculated. 

Comparisons  of  zone-size  and  fracture-stress  calculations  with  experi¬ 
ments  show  that  predictions  of  the  modified  DM  model  can  be  reasonably  ac¬ 
curate.  Furthermore,  it  is  a  powerful  tool  in  the  hands  of  the  alloy  devel¬ 
oper,  since  it  shows  quantitatively  how  much  strength  level,  ductility,  and 
the  work-hardening  rate  contribute  to  fracture  toughness.  Preliminary  cal¬ 
culations  offer  promise  for  describing  the  propagation  of  ductile  cracks  and 
the  transition  from  ductile- to-brittle  fracture.  These  calculations  are 
limited  by  the  lack  of  good  stress-strain  data  at  strain  rates  approximating 
those  at  the  tip  of  a  fast-running  crack. 
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APPENDIX  A:  NOTATION 


Crack  plus  plastic-zone  length 

Parameters  in  the  stress-strain 
equations 

Crack  length 

Point  in  plastic  zone  where 
step  change  in  internal  tension 
is  applied 

Plastic-zone  width 

Young's  modulus 

Fracture  toughness 

Strain-rate  sensitivity, 
d  log  o/d  log  e 

Strain-hardening  exponent 

Reduction  in  area 

Internal  tension 

External  applied  stress 

Crack-extension  stress 

Specimen  thickness 

Minimum  specimen  thickness  at 
crack  tip 

=  Crack  velocity 

=  Trial  crack  velocity 

=  Displacement 

=  Crack-tip  displacement 

=  Distance  from  center  of  crack 
in  direction  of  crack 

=  Discrete  value  of  x  used  in 
computation 


y  =  Coordinate  axis  normal  to 

crack  plane 

Y  =  Yield  stress 

a  =  arc  cosh  x/a  (x  ^  a) 

j3  =  arc  cos  c/a 

e  =  True  strain 

7*  =  True  strain  at  fracture 

j  =  Plastic-strain  rate 

e  =  True  strain  at  crack  tip 

c 

7  =  L'uders'  strain 

L 

y  =  Poisson's  ratio 
0  =  arc  cos  x/a  (x  £  a) 


eJ 

= 

arc  cos  Xj/a 

(xj  £ 

a) 

p 

= 

Plastic-zone 

stress) 

size 

(plane 

Ph 

= 

Plastic-zone 

strain) 

size 

(plane 

a 

= 

Stress 

0 

= 

True  stress 

T 

= 

Time 
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APPENDIX  B:  THE  DISPLACEMENT- STRAIN  RELATION 


A  simplified  model  of  the  plastic-zone  strain  distribution  is  shown  in 
Figure  B-l.  Figure  B-la  sh  s  the  cross  section  of  a  plate,  of  thickness  t, 
at  the  distance  c  ^  x  £  p;  the  line  AA  identifies  the  plane  of  the  crack. 

The  deforming  region  is  arbitrarily  confined  to  a  height,  d.  It  is  also 
assumed  that  strains  in  the  x-direction  can  be  neglected.  Under  stress, 
this  region  extends  to  a  length  (d  +  2v)  at  constant  volume,  and  the  plate 
consequently  necks  as  shown  schematically  in  Figure  B-lb.  Necking  tends  to 
localize  deformation,  but  the  neck  cannot  become  too  small  or  plastic  deforma¬ 
tion  will  be  constrained.  Although,  to  account  for  constant  volume  deforma¬ 
tion  the  plastic-zone  width  grows  larger  than  d,  for  the  majority  of  cases 
treated  here,  this  growth  is  small,  d  »  2v.  Thus  the  width  of  the  plastic 
zone  remains  roughly  constant.  The  neck  is  approximated  by  the  double 
trapezoid  of  Figure  B-lc,  Since  the  deformation  involves  no  volume  change, 
the  area  of  the  trapezoid,  d*t,  is  constant,  and  this  fixes  the  minimum 
section  t' , 


The  maximum  true-strain,  e  ,  is  the  strain  corresponding  to  the  minimum  sec¬ 
tion, 


The  tension,  S,  supported  at  x  and  o,  the  true-stress  acting  at  the  minimum 
section  can  also  be  related, 

s  "  °  (t")  =?  (di-^) =  ®  exp(-?)  •  (B_3) 

In  this  way,  values  of  ?  o,  and  S  can  be  calculated  if  v,  d ,  a  =  f(e)  are 
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